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1.0 INTRODUCTION: 


This Semiannual Status Report summarizes the work performed under Grant 
NAC5-458 entitled "Investigation of Dynamic Noise Affecting Geodynamics In- 
strumentation in a Tethered Subsatellite." 


During the reporting period, SAO has 

- Written an invited paper to be polished in the forthcoming Special 
Issue on Geodynamics of w .he IEEE Transactions on Geoscience and Remote 
Sensin g . 

- Continued modeling of the atmospherically induced dynamic noise in the 
subsatellite through the modification and use of the SKYHOOK program; 
interaction wiuh an elementary form of random atmosphere was simulated. 

- Developed a method for stabilizing the subsatellite against the rota- 
tional effects of atmospheric perturbations and analyzed a simplified 
model . 

- Performed a variety of analytic studies of tether dynamics aimed at 
elucidating dynamic noise processes and at continuing the rotational 
stabilization analysis. 

- Discovered a novel mechanism for coupling longitudinal and latitudinal 

oscillations of the tether: a longitudinal (approximately vertical) 

oscillation causes an oscillation in the subsatellite's altitude; this, 
in turn, causes an oscillation in the drag experienced, which is in the 
horizontal (latitudinal) direction. 

- Begun study of random vibration analysis for modeling the TSS under 
atmospheric perturbation . 
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2.0 INVITED REVIEW PAPER 

A paper entitled "Gravity Gradiometry from the Tethered Satellite System", 
by G. E. Gullahorn, E. Fuligni and M. D. Grossl, was written and submitted for 
the Special Issue on Geodynamics of IEEE Transactions on. Geoscience and Remote 
Sensing. This paper has been accepted and is scheduled for publication in the 
July, 1985 issue. It has been modified substantially from the draft version 
included as an appendix to the Semiannual Report. 

Our survey of the open (published) literature showed virtually nothing on 
the topic of TSS gravity gradiometry; even among less readily available sources 
such as contract and grant reports or limited distribution summaries of meet- 
ings, nothing at an introductory level was found. Thus, in keeping with the 
spirit of a "Special Issue" intended for a broad audience, we decided to write a 
largely tutorial/review paper rather than simply report one or two current 
results comprehensible only to a few other specialists. We did, however, also 
summarize some SAO work in the field. 

Although researching and writing this paper was a substantial effort, we 
believe that the product was valuable. It will publicize the concepts and 
advantages of TSS gradiometry (and of space-born gradiometry platforms in gen- 
eral) , and provide a ready entr6e for the fledgling researcher as well as the 


Interested scientist from a related field. 
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2.1 The Significance Of Platform Rotation 

The discipline of producing such a review was alao of value to the authors. 
The process of TSS gradiometry, and the role of dynamic noise calculations, was 
placed in context. The important role of rotation of the instrument platform 
(subsatel lite) was examined and emphasized. Excerpting some paragraphs from the 
paper : 

The gravity gradient is a tensor whose components in Cartesian coordinates 

are defined as 


Tij = dgi/dx 3 = - d 2 V/dxidx] (2.1.1) 

The gravity gradient is expressed in Eotvos units, 1 E = 10'* gal/cm = 10'* 
sec'*. Note that because the gravity vector arises from a potential, Tij = Tjt; 
and in free space Poisson's equation will relate the three diagonal components, 
Txx * lyy + r„ = 0. Thus, of the nine components, only five are truly indepen- 
dent . 


Measurements made with an Instrument linearly accelerated relative to a 
rest frame will not be affected due to the phenomenon of common mode rejection 
discussed below. However, if the instrument is subject to a rotation with 
angular velocity £3( t) , the measured gradient will be 


paaaa _ pr«»t -f T + T J 


( 2 . 1 . 2 ) 


where r r#at is the gradient of (2.1.1) , T(t) is the matrix 
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T(t) 


0 

-u>, 

w 2 

Wj 

0 

-III! 

-w 2 

Wl 

0 


(2.1.3) 


and T J means matrix multiplication in the common sense. Note that the term 
dependent on the angular acceleration ic anti-symmetric and can be immediately 
distinguished from the inertial frame gradient if the full tensor gradient is 
measured accurately enough; the term dependent on the angular velocity is sym- 
metric, and cannot be removed from the observations without knowledge of the 
instrument's rotation. 


If the gradiometer is not attached to a completely motionless platform, the 
measured gradient may contain errors. A linear acceleration, with no rotation, 
is relatively innocuous due to the phenomenon of "common mode rejection": the 
two test masses will be subject to the same non-gravitational acceleration, so 
the difference, the measured gradient, will not be affected. Rotation of the 
instrument, however, will Innately alter the sensed gradient as shown by 
(2.1.2). If the rotation has angular velocity w(t), ignoring for simplicity the 
vector nature, then the T J term of (2.1.2) will have magnitude approximately u 2 . 
Measurements at the goal level of 3 x 10 4 E = 3 x 10" 11 sec' 1 will be degraded 
if u > 5.5 x 10' 7 sec 1 - 0.1 deg/hr. Stability at this level should be 
readily achievable, and if not, such rotational velocities can be easily 
measured and their effect removed via (2.1.2). The rotational acceleration 
term, T, however, will interfere with measurements if u > 3 x 10' 1 * see'*; to 
subtract the effects would require, over a one second integration, knowledge 
of u to about 3 x 10* 1 * sec' 1 = 6 x 10' B deg/hr, well beyond current technology. 
If eight of the nine components of T are measured, the rotational acceleration 
effects can be removed by their anti - symmetry . However, even a very small u 
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will generate large spurious gradients, and an enormous dynamic range require- 
ment, so it will still be very important to isolate the gradiometer from rota- 
tional acceleration. Alternately, one can measure the diagonal terms only, 
which are not affecred by the rotational acceleration; but then one must be 
careful to avoid contamination by the cross terms, due, for instance, to 
misalignment of the accelerometer pair. 

Another source of error lies in the vertical 3000 E background due to the 
overall Earth mass. Any slight misalignment of the spacecraft reference frame 
will project this on the otl sr axes, e.g. 1/2 arc sec will lead to 10' J E. 
This is not critical in itself, but any fluctuation in alignment at the arc 
second level could cause serious problems . 

3.0 NUMERICAL SIMULATIONS 


The numerical simulations being done at SAO use a modified version of the 
SKYHOOK program, which models the TSS as set of discrete masses and connecting 
(massless) tether segments. As described in the Final Report to the previous 
contract, NAG5-325, we have added facility for superimposing spatial fluctua- 
tions on the standard atmospheric density routine and outputting the resultant 
accelerations of the subsatellite, together with ambient densities at each mass, 
for later analysis. The acceleration components found to be most elucidating 
are those tangent to the tether at its attachment to the subsatellite and or- 
thogonal. Expressed this way, the effects of tether longitudinal vibrations are 
separated from the direct effects of atmospheric drag variations. Restrictions 
are that the subsatellite is considered a point mass (attitude dynamics and 
aerodynamic effects are not included) and that computational effort limits the 
number of masses to about 10 for routine use, hence limiting the frequency 
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response of the system. All cases considered have beer, limited to a tether 
deployed in the orbital pl.ane, though this restriction is not inherent in the 
program . 

3.1 Software Modifications 

3.1.1 Graphic Capability Enhr ncement - 

One simple, yet important, enhancement of SAO capabilities was the adapta- 
tion of graphics tools allowing informative and compact display of the accelera- 
tion and spectral outputs. Compare the figures in this report with the hand 
drawn Figure 7-1 or the printer plots in Appendix D of the NAG5-325 Final 
Report . 


3.1.2 Modification Of Underlying (Mean) Atmospheric Density Routine - 

We resolved a subtle problem in the standard SKYHOOK atmospheric density 
routine which had been causing an occasional sharp jump in the orthogonal accel- 
eration component. The difficulty arose from an interaction of two factors: 
First, to eliminate residual deployment effects we have been forced to make two 
SKYHOOK runs for each case considered, a "reference" run with no atmospheric 
perturbation and a run with the perturbation we wish to ..tudy. The resulting 
accelerations are differenced before plotting, computing spectra or other analy- 
sis. Second, the atmosphere routine used simple linear interpolation in a 
table, resulting in a p (h) curve with only piecewise continuous derivative. 

Because the two runs are different, the subsatellite will be at slightly 
different altitudes in each and thus subject to slightly different drag. If 
both subsatellites are in the same tabular Interval, thla difference will be 
constant or vary only very slowly, causing minimal effect. But when one sub- 
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satellite croatas a tabular altitude, and thus J.s subject to a different density 
as function of altitude, this difference in drag takes a sharp jump. 

To remedy this problem, we created a smoother atmospheric density routine. 
We first attempted to simply interpolate log(p), but ended up fitting a moderate 
order polynomial to the log(p) data (with independant variable 
^altitude- 10 Okm ) at each tabulated tempo) ature and Interpolating in tempera- 
tur« . This produced the desired smoothness with no noticeable deterioration in 
execution time. Our scheme of differencing the two runs may still have some 
effect at very low frequencies, as the relative altitudes of the subsatellites 
change gradually over the run. 

Tracking down this effect of altitude on drag did point up one interesting 
mechanism : 


Any variation in the subsatellites altitude due to tangont 
accelerations produced by the tether (approximately vertical) 
will cause a variation in the ambient atmospheric density, 
hence the drag experienced, hence the orthogonal (approxi- 
mately horizontal) acceleration. 


3.1.3 Multiple Smooth Perturbations - 

The version of SKYHOOK created during the previous contract had allowed 
only a single density perturbation, with a sharp cutoff between the enhanced (or 
diminished) region and unperturbed density; the region is elliptical in cross 
section and extends indefinitely to either side of the orbital plane. Apart 
from the obvious restriction to a single perturbed region, the sharp boundary 
results in very long integration time due to the peculiarities of the Gear 
integrator used to solve the differential equations of motion. 
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W« hav* extended th* allowed perturbations to include up to two hundred 
(possibly overlapping) regions each with a smooth cutoff to enhance computa- 
tional efficiency. The parameters of each region are: a vertical radius, r v ; 

a horizontal radius, r h ; altitude of center, h c ; a distance "along orbii " from 
the start of the simulation, a e ; and an enhancement factor, e. Again, since we 
are concerned with deployment in or near the orbital plane, the density per- 
turbations extend indefinitely perpendicular to the plane. First, we define a 
scaled radius 

r « y/ ( (h-h«) /r v ) 1 ♦ ((a-a c )/r h ) J 

where h is the altitude and a the "along orbit” distance. Then ihe density 
perturbation for a single region, relative to the "base" unperturbed density, is 

P/Pb.— = 1 ♦ e/(l*rV. 

The region so defined has the same total mass (in a two dimensional sense) as 
the previous sharply bounded region with the aame vertical and horizontal radii 
and enhancement factor. When there is more than one perturbed region, th* right 
hand side of this equation is computed for each perturbation and they are all 
multiplied together to form the total perturbation. 
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9.2 Simulations Of Single Psrturbations 

Ws havs analyzed four casen: smooth and sharp bounded regions; impinging 
on ths subsatellite at 120 km and on the tether at 170 km altitude. All regions 
were circular with 20 km (effective) radius., and 20% enhancement. 

A typical plot of the residual tangent and orthogonal accelerations experi- 
enced after encounter with the perturbation in given In Figure 3-1; the other 
cases were similar. Note tho smoother nature, and much smaller magnitude, of 
the orthogonal acceleration. 

The spectra of the accelerations produced in the various cases are also 
surprisingly similar. Figure 3-2 shows those for two cases: a smooth perturba- 
tion Impacting directly on the subsatellite, and a sharply bounded perturbation 
impacting halfway up the tether. The sharp peaks corresponding to the spring- 
mass mode and the expected eight longitudinal vibration modes are clearly 
visible in the tangent acceleration. These unperturbed modes are mimicked by 
the orthogonal acceleration, with the addition of a mode at very low frequency, 
possibly due to latitudinal vibration modes. The spectra are largely similar, 
though the peaks in the case with tether Impact are smaller in magnitude; note 
that the background levels are the same in both cases . To compare in more 
detail we tabulated the frequencies and magnitudes of the spectral peaks from a 
printout. The frequencies are all very close, to within ±0.001 Hz. The mag- 
nitudes are plotted In Figure 3-3. The disturbances with smooth and sharp 
cutoffs, at subsatell kte altitude, are very similar in both magnitude and rela- 
tive strengths of the pesks. The disturbance at mii-tether altitude also ex - 
elves the modes in roughly the sam«. ratios, although with substantially less 


magnJ tude . 
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Fig . 3-1 . Acceleration* , in gals, produced by a single atmospheric density 
enhancement impacting the subsatellite. Smoothly cutoff region, 20 ka radius, 
20% enhancement. Acceleration as a function of tiae is plotted, the component 
tangent to the tether on top and orthogonal on bottoa. Impact with the enhanced 
region occurred before the data plotted, although soae of the final direct ef- 
fect is visible in the orthogonal component . 
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Fig. 3-2 . Logai lthm of the acceleration spectral magnitude plotted as a func- 
tion of frequency (Hz). In each plot, the top spectrum is for the tangent 
component, the bottom for the orthogonal. The top plot is for the case illus- 
trated in Figure 3-1, a single smoothly falling off enhancement impacting di- 
rects* on the subsatellite. The bottom plot ie for an equivalent region impact- 
ing on the tether at ' ts center, sharply cut-off to avoid any direct impact on 
the subsatell « . 
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SPECTRAL PEAKS COMPARE SMOOTH AND SHARP BLOBS AT 120 KM ALTITUDE 




Fig. 3-3 . Comparison of the magnitudes of spectral peaks. The peak magnitudes 
are on logarithmic scales. The top graph compares equivalent regions with sharp 
and smooth cutoffs, both impacting directly on the subsatellite. The bottom 
compares identical enhancements impacting on the subsatellite ar. 120 km altitude 
and on the tether only at 170 km altitude. In each graph, the top pair of lines 
is for the tangent component, the bottom for the orthogonal component. 
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3.3 Modelling A Random Atmosphere 


Simulations of encounters with individual "bumps" in th« atmospheric den- 
sity serve to elucidate the impulse response of the TSS . However, it would also 
be instructive to simulate the TSS passing through a realistic randomly varying 
atmosphere . 


From Gross, Reber and Huang (1984) and Gross and huang (1985) , we know that 
the power spectral density of variations at 250 km is roughly a power law with 
exponent -3 at low spatial frequencies and -3 to -4 at higher frequencies. This 
appears to be the best information available, no similar data existing at the 
lower altitudes of primary interest to us. The best that can be done, in the 
absence of some physical theory, is to simply assume the same per-cent variation 
relative to mean ambient density in the regions of interest. 


Creating a random atmosphere by superposing a set of the perturbed regions 
discussed above was attractive for several reasons: 

- By analogy to individule "turbules" in a turbulent medium. 

- The individual elements can be implemented and studied before combining 
them. The superposition is then a relatively simple programming task. 

- The Gear integrator used in SKYHOOK demands smooth data for efficient 
operation. Interpolating on a grid, for instance, would cause much step 
size and integration order adjustment crossing each mesn line. 

- The amount of data needed to represent the atmosphere at a one second time 

resolution in the horizontal direction and an equivalent distance (say 8 
km) vertically is large: roughly 50,000 points for a one hour run. 

- Because the Gear integrator, ir adjusting step size, may look at intermedi- 
ate time values after taking one large time step, the density must be 
effectively created in toto rather than step by step as the simulation 
proceeds. Having the density at the same place and time change in value 
depending upon the series of time steps used would be likely to wreak havoc 
with the integrator's performance. 


We have started an investigation of how to create a set of perturbed regions for 
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use by the current version of SKYHOOK which will have a given spectrum (as seen 
by the subsatellite traversing the atmosphere). This method apoears, in light 
of these efforts, not to be capable of producing the desired atmosphere model 
(an individual realization of the random atmosphere) . 

Two methods appear to show promise, though time was not available to pursue 
them: First, as an ad hoc approximation assume that the masses in SKYHOOK are 
at constant altitude. For each mass, create a time series with the desired 
resolution, interpolate this with a high order spline, and use these separate 
series. This would still require a substantial amount of data, but would be 
less strongly interconnected and require no interpolation in the vertical direc- 
tion. Although acceptable, perhaps, for SKYHOOK with it's limited number of 
masses, this method could not be extended to simulations requiring greater 
vertical resolution. Second, the "turning lines method" (see Bras and 
Rodriguez- I turbe , 1985) allows reconstruction of a two or three dimensional ran- 
dom field from distributions along a set of lines through a common center; the 
number of lines suggested is small, 4 to 16. This method is for producing a 
distribution on a grid, but may be adaptable to our need for densities at unpre- 
dictable points. 

In the absence of a good theory, we have created a set of 200 regions with 
random centers (altitudes 120 to 220 km, random distance along track), random 
radii (up to 30 km vertical, 50 km horizontal) and random enhancement factors 
(in the range -0.5 to +0.5); all distributions being uniform. The induced 
accelerations are shown In Figure 3-4, and the spectra in Figure 3-5. The 
tangent acceleration is roughly similar to that produced by single regions, and 
in the spectrum we see the same peaks as before although more weakly (note the 
higher baseline) The same modes are being excited, but with more "noise" 
generated by the random forcing function. The orthogonal acceleration, however, 
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Fig- 3-4. Accelerations produced by a model atmosphere with 200 randomly 
scattered enhancement regions . The orthogonal component (top graph) appears to 
be responding directly to the ambient atmospheric perturbation, whilti the tan- 
gent component (bottom graph) is more similar to the results from a single 
region model (Fig. 3-1) . 





OHGWAL PAGE iS 
OF. POOR QUALITY 


Page 18 




0 06 .1 16 2 26 3 .36 .4 46 6 

FREQUENCY (Hi) 


Fig. 3-5. Spectra of the time series in Figure 4. Note that the tangent compo- 
nent spectrum shows distinct peaks corresponding to exalted vibration modes, 
while the orthogonal component spectrum is largely featureless. 
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Is totally dominated by the direct effect of the drag variations, and its 
spectrum is featureless. For reference, the ambient density at the subsatellite 
is shown in Figure 6, along with its spectrum. These are similar in character 
to the orthogonal acceleration, as expected if we are seeing largely direct drag 
effects in the latter. 

4.0 ANALYTIC INVESTIGATIONS 

SAO has performed several analytic Investigations of tether and TSS system 
dynamics. These fall in two categories: those aimed at overall understanding 
of tether dynamics and it's implications for dynamic noise, and those specifi- 
cally aimed at developing noise abatement methods. We treat the latter first. 

4.1 Simple Model For Eliminating Rotational Perturbations 

As pointed out in Section 2.1, rotational perturbations of the instrument 
platform are much more insidious than simple translational accelerations. The 
SKYHOOK simulations described in Section 3 do not Include subsatellite rotation, 
but it is easy to see that intersection with a density enhancement must cause 
the subsatellite to rotate: the center of force will be near to the center of 
the subsatellite, rot through the attachment point to the tether, and hence any 
differential acceleration of the subsatellite and tether will cause a torque. 

To alleviate this torque, we adopt the idea of a ballast mass: a second 
mass attached by a short piece of tether (or rigid rod) to the instrumented 
subsatellite and with an area/mass ratio chosen so that drag variations will not 
excite rapid oscillations. We assume that this suspended subsystem is small 
enough thrt all points experience the same atmospheric density at any one time, 
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Fig. 3-6 . The ambient atmospheric density experienced by the subsatellite in 
the simulation of Figures 3-4 and 3-5, and its spectrum. Note the similarity to 
the orthogonal component's graphs in Figures 3-4 and 3-5. 
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and that th* gravity gradient strength (which w* shall simply call "g") Is 
uniform. 


A full analysis of such a system, some attempts at which are given In the 
next section, promises to be exceedingly complex, probably being Incapable of 
analytic solution. Here, we present a very simplified model problem which, 
however, demonstrates the feasibility of the ballast mass Idea and Illustrates 
the basic form which a complete engineering analysis would probably follow. The 
ballast mass concept, and our simple model, are illustrated in Figure 4-1. We 
model the system as a pair of coupled pendulums suspended from an "infinite 
mass" Shuttle In the orbital plane, and assume small deviations so that the 
system becomes linear. We take as the reference configuration simple hanging 
equilibrium without the complication of the equilibrium displacement due to the 
mean drag, which should be adequate for the nearly vertical configurations ex- 
pected. The tether segments are assumed massless and inextensible , and the 
masses are point masses; the rotation will be apparent in the motion of the 
angle of the segment joining the two masses. The equilibrium tensions in the 
two segments , which we may assume to be constant , are 


Ti = g(mj. + m 2 ) 


(4.1.1) 


T a = g m 2 

The system state is fully described by the two angles flj. and 0 2 . We shall 
explicitly balance forces in the x (horizontal) direction, and express the 
result in 0 2 and 0 2 . The homogeneous part of the equations is also easily 
derived via a Lagrangian, leading to the same result. We have: 

m 2 Xi = — Ti sin (0 2 ) + T 2 sin (0 2 ) + d 2 


(4.1.2) 


miXi = -T a sin(0 2 ) + d 2 


To evaluate the left hand sides, use 


xi = L sin(0i) = L 0i 


(4.1.3) 
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Fig. 4-1 . Simplified model of the ballast mass concept for alleviating rota- 
tional perturbations due to atmospheric density variations . The system is 
modelled as two coupled pendulums, with inextensible massless tether segments 
and point masses . 
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Xa = xi + b sin (f 2 ) = Lf x + b fa 

and similarly expand tha sin() functions on the right hand sides. Use the 
dimensionless parameters 

e = b/L 

a = dx/da = Ax/Aa (4.1.4) 

P - ®a/mx 


where Ax and A; are the cross sectional areas of the two masses; e will be a 
small parameter, while a and P may be manipulated to achieve desired behavior 
and are of order unity. Also write d = da and m = n>x , and after some manipula- 
tion obtain: 



1 o' 



L 

1 « 

. *». 

+ 9 


1 +P ~P 

0 1 



V 


a 




. V*. 


d(t) 

m 


(4.1.5) 


where we explicitly show the dependence of drag on time, d(t). In more concise 
matrix notation, 

Al + B# = a d(t) (4.1.6) 


To solve this system of equations we shall assume that it is initially at rest 
and apply the Laplace transform 

£[f(t)] = F(r) =^e~f(t)dt (4.1.7) 


Using the properties 

£[*' (t)] = - f (0) + rF(r) 
r[f"(t)] - - f ' (0) - r f (0) + r J F (r) 
along with the initial conditions, and letting 

D(r) = £[d(t)] 

0(r) a [6i • 0,]* = £[#(t)] 

we find 


(4.1.0) 


(4.1.9) 


r’A0 + B8 = aD (r ) 


(4.1.10) 
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This is an algabrale aquation for 0(r) , and can b* aolved to giv* 


Ox 


1 

A 


Lr , «'»g qfi 

-Lr* Lr a +q{l+P) 



a 


Wfi. 


Lirl 


where 


A s (Lr* + g(l + 0)) (Lr J « + g) + qfiLr* 

Wa a ra spacifically intaraatad in 0 2 : 

e,(r) = < Lf ‘ + 9(1 +<»))/? - «Lr- D(r) 

A m 


(4.1.11) 


(4.1.12) 


(4.1.13) 


To progreaa further, v* serum* a simple impulse, d(t) = d« 6 (t-T) . 17, an 

D(r) = d#*‘ Tr . For th* moment, let d« = m = 1 ; if needed, linearity of the 

Laplace transform would lat ua factor tham in after tha answer is obtained. 
Also, we can let T = 0 , since the factor *' ,r in the Laplace transform simply 

shifts th* time-domain function by T: £‘ l [*‘ TT F(r)] = f (t-T) . Thus, w* may 

take D(r)/m = 1 in (4.1.13). 


We now need the inverse Laplace transform of th* (known) 0 2 (r). From here 
on, MM shall aaauma t « 1 , that if h ^ L- Then we factor the denominator to 
lowest order in < : 


« L 1 (r 1 + * a ) (r* -f f») 


(4.1.14) 


where 


,i = gt 1 +P ) = 3 1 1 +J) 

Lc b 


(4.1.15) 


Aa we shall see below, f represents a alow mode of overall oscillation, while a 


gives a fast mode of subsatellite oscillation. Noting that 

1 _ 

(r» + *>) (r* + f») " f 1 

and that f*-e* ■ -** to lowest order in e, and letting 


b [ Trb - ] <4116) 


h l . 

W-aft) 


(4.1.17) 


we find after considerable algebra 

e,(r) 


- [wrffc] <4 l l8 > 
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Putting th* factor d^O-L/m back into th* answer, and using £' x [ 1 / ( r 2 *a 2 ) ] * 

(1/a) sin (at) , w* obtain 

»*«•> - - [ £ ] [ ] { ( * 4 ^ ) ■ ln! * t) ( * 4 ^ ) * ln(,t) } < 41l9 > 

W* hav* two parameters, a and 0, at our disposal and ws us* them to elimi- 
nate the fast mode corresponding to sin (*t) by setting h 2 -* 2 = 0 , which leads to 
a0 “ 1-c , or 

= ( 1-0 U (4.1.20) 

to first order in c . 

This answer is in some sense obvious, to zero'th order in c: for the 

masses to accelerate equally, neglecting the suspending tether, their mass/area 
ratios must be equal. What we have done here is to derive a first order approx- 
imation and, more importantly, illustrated techniques that might be used in a 
more complicated and realistic calculation. The fact that we have been able to 
eliminate the fast mod* in this first order calculation, as well as the zero'th 
order one, also is promising for success in eliminating fast rotational modes in 
more complete models. 

Before leaving this very simple model , we should consider one last ques- 
tion: What is L? Most simply, if the assumption cf a massless, inextensibl* 

tether were precise, L would simply be the distance from the subsatellite to th* 
shuttle. The real tether, however, will hav* an influence and an L useful for 
back-of- the- envelope calculations might be taken as, say, the length of tether 
whose mass equals the subsatellite (not substantially different in the case of a 
100 km tether with density 8 kg/km and a 500 kg suksatellito) , or the wavelength 
of transverse waves having th* same period as the subset-ellite pondular oscilla- 
tion (which is much smaller: for a subsatellite with 0.2 Hz oscillations and a 

tether with spectrum O.OOln Hz, this gives soda n=200, or wavelength about 1/2 
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ka) . Providing a confidant anivar to this quaation will dapand on axparianc* 
froa a amp la calculations including a continuum tet K *r; th* baginnings of such 
calculations ara discussad in tha naxt section. 

4.2 Extansions To Tha Siapla Modal 

In tha abova saction va have domonstratad tha faasibility of using a bal- 
last aass to halp aliminata rotational parturbatlona of tha subsatallita instru- 
ment platform by atmospheric density variations. A great deal remains to ba 
dona before such abatement methods can ba considered fully designed. 

Still within th* context of non-continuum tethers, on* could add to the 
modal such factors as tether elasticity and damping, rotational dynamics of th* 
subsatallita and ballast, and more complicated arrays of ballast masses (e.g., 
ons below and on* abova th* instrument platform) and specially designed connect- 
ing elements (such as highly damped cables) . Such compi.i itions are likely to 
soon reach the limits of analytic solubility and require numerical solution and 
various numerical parameter optimization methods. 

Adding the continuum tether adds another whole level of complexity, requir- 
ing solutions to mixed partial and ordinary differential equations, but must be 
done before th* problem can be considered fully explored. We have achieved a 
partial solution of s simple case: one point mass on the eno of a tether with 
no dampir^ and no body forces (drag, gravity) , hanging from an infinite mass 
shuttle. The suspended aass experiences a variable drag d(t). 

We shall not give our calculations in detail. The oscillations in th* 
three orthogonal directions are decoupled (for small displacements) , and only 
the horizontal, along-orbit direction is effected by d(t). The equations for 
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this component are: 

MVtt = (T.f/L) V*b 0 < 8 < t 

v (0 , t) = 0 
v(/,t) = y\ (t) 

m fj (t) = - (T.f/L) v.(/,t) + d(t) 

where t is the tether's natural length, L its stretched length in equilibrium, 
T. the equilibrium tension, »j (t) the displacement of the end mass, and the rest 
of the notation follows previous discussion. Laplace transforming this system 
with respect to time generates an ordinary differential equation with r as a 
parameter : 

r*V(s,r) =- c* V.. (s,r) 0 < s < l 

V(0,r) = 0 
V(f,r) = E (r) 

r* E(r) = -b* V.(f.r) t D(r) 

where V and E are Laplace transforms of v and rf , D is the transform of d(t)/m, 
and b and c are constants. This system is readily soluble to give 


V(s , r) 




D(r)/r« 

slnh(rf/c) + (b*/rc) cosh(rf/c) 


sinh (rs/c) 


E(r) 


= 2112 

r* 


sinh ( rf/c ) 


sinh(rf/c) + (b*/rc) cosh(rf/c) 


We are primarily interested in the motion of the mass and thus in the inverse 
transform of E(r). Unfortunately, even for simple forms of D(r) this does not 
seem to be tabulated nor to be modifiable to a recognizable form. Further 
investigation may provide an analytic solution, but fully realistic calculations 
with systems of suspended, finite extent masses will likely require numerical 
approaches. A modal (separation of variables) analysis should also be made to 
attempt to find the response to a drag impulse. 
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4.3 Cooperation With Prof. Silvio Bergamaschi 

In December 1984 Prof. Silvio Bergamaschi of the Institute of Applied 
Mechanics in Padova, Italy visited SAO. Though NAG5-458 provided no direct 
funds for his visit, we did spend several days discussing tether dynamics, par- 
ticularly dynamic noise, with Prof. Bergamaschi. Prof. Bergamaschi is develop- 
ing a model of the TSS including a full continuum tether with end masses, and 
gravity gradient force on the tether itself. Modal anal/sis is being applied to 
this model. We assisted him with seme numerical mode calculations, and dis- 
cussed some non-intuitive physicul implications of the results. We hope to 
continue the cooperation in developing this model, which will provide both a 
valuable check on SKYHOOK numerical computations, and results not obtainable 
with SKYHOOK. In the future we hope to cooperate with Prof. Bergamaschi on 
measurement and analysis of the TSS dynamic environment during the TSS Demon- 
stration flights. 

4.4 Effects Of Tether Internal Damping 

Our simulations and analytical calculations have ignored the effect of en- 
ergy dissipation within the tether. The simplest form of such dissipation would 
be Internal viscous damping. This has been partially because of the added 
complications this would introduce and partially because of uncertain knowledge 
of tether properties. We have, however, made some simple theoretical calcula- 
tions which we present here and which indicate that tether viscous damping may 
be quantitatively significant though probably not qualitatively crucial to un- 
derstanding of dynamic noise. (Some of the work in this section has been per- 
formed in conjunction with NASA Contract NAS8-36160, and a portion of this sec- 
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tlon has been excerpted from Quarterly Reports 2 and 3 of that contract.) 

We begin at the very beginning, and derive the equations governing a tether 
with internal viscous damping. Though certainly not new, this material does not 
seem to be readily available (typical sources give equations for cables subject 
to external damping, such as an ambient fluid) . First we consider a tether 
se gment . for the moment considered massless; later we shall derive the partial 
differential equation for a continuum tether. We wish to derive the dependance 
of the damping constant for a piece of tether material on the length of the 
piece, under some very general assumptions. Specifically, we assume that the 
material is viscous but not plastic . That is, there will be resistance to 
motion, but no permanent change in the equilibrium state. We assume a simple 
Hooke's law form for the elastic portion of the stress, and the viscous damping 
will add a term proportional to velocity. 

Consider a simple physical system: a length of tether fastened to a wall 
at one end and a mass M at the other. Neglect the tether mass. Let x be the 
extension of the segment past natural length. Then the equation of motion is 

Mx = -Kx -Bx (4.4.1) 

where K is the Hooke's law spring constant and B is some damping constant. Now 
suppose we cut the tether in its center and place a small mass m between the 
pieces; denote the extension of the segment between the wall and mass m by y. 
Then the equations of motion are twofold: 
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Mx = -k(x-y) -b(x-y) 

my = - ky -by (4.4.2) 

+ k (x - y) + b (x — y) 

If we let £ = (x-y) , a bit of manipulation arrives at the equation 

Mx = - %kx — %bx - %m(x-£) (4.4.3) 

Let m -* 0 and compare the result to equation (1). If the accelerations in 
the third term of (3) remain finite, which we shall simply assume since this is 
not meant to be a rigorous proof, agreement between these two results for the 
motion of mass M requires that k = 2K, b = 2B. Generalizing the argument, we 
have that 

k ex i bay (4.4.4) 

where i is the (natural) length of the tether segment. Tor the spring constant 
k, the proportionality constant is simply EA, Young's modulus times the cross 
sectional area. For the damping constant b, we shall follow a remark in Bodley 
and Park (1983; p. 54) and denote the proportionality constant by G*. (Note 
that EA and Cy apply to a tether of given cross section and manufacture. It 
seems raasonable that the energy dissipation, hence the resistive force, will be 
proportional also to the cross sectional area of the tether; we may expect that 
c v = Cv/A is a material constant analogous to E.) 

Now let us derive the partial differential equation for longitudinal motion 
of a tether with viscous damping, considered as a continuum. We suppose that 
the tether remains in a straight line, and is sufficiently extended in its 
equilibrium state (e.g., by loading) so that it does not go slack for small 
perturbations. Let u(s,t) be the displacement from equilibrium of a mass ele- 
ment which would be at coordinate s in equilibrium; thus the full position will 
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be x(s,t) = s ♦ u(s,t). The tether has properties AE and Cy as above, and mass 
per unit length p. 


Apply the lumped mass or ball-and-spring approximation to the continuum 
tether, with masses separated by As = h. The mass of each particle is then ph, 
and denote it's displacement from equilibrium by Ui (t) . The acceleration of a 
particular mass will then be given by the sum of the spring forces and damping 
of the two tether segments on each side; these will each have k = AE/h and b = 
Cy/h. 


MhU t = ♦ (AE/h) (Ufi-Ui) + (C v /h) (U w -UO 

- (AE/h) (Ut-Ut-0 - (Cv/h) (Ui-Ui-i) 


(4.4.5) 



We are now able to examine the effects of damping. Consider a tether of 
length L, with fixed ends as a "standard" case. With a = AE/fi and b = Cy/p, 
the PDE becomes u t t = (a ♦ bd/dt) u,. . Standard separation of variables with 
u(s,t) = S(s)T(t) yields modes 


S„(s) = Sin (XnS) 
T„(t) = exp (a (*„) t) 


(4.4.8) 


where 


a(K) 


k„ = n r/L 


(1/2) [ -K'b ± 



(4.4.9) 


The wavelength of a mode associated with «c is just W = 2 w/k. 
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Note that for K , corresponding to the longest wavelength modes, the square 
root in (4.4.9) in Imaginary and the modes will be oscillatory. For non-zero 
damping, though, the higher modes will be overdamped with no oscillatory behav- 
ior. This occurs when the square root becomes real, at 


Worn E » C» / y/ fiAE 
ncait — 2L/Wan 


(4.4.10) 


Note that this critical wavelength for damping depends only on the tether 
properties and not on the length of tether involved; the total length deter- 
mines only how many oscillatory modes will actually exist. Of related interest 
is the maximum achievable frequency of an oscillatory mode. The frequency, that 
is the imaginary part of o(x), is zero for #c=0 , rises to a maximum value, and 
declines to zero at * = 2ir/Wat IT ; above this, the modes are not oscillatory. 
This maximum value is 

= AE/Cv (4.4.11) 


To evaluate these quantities we need estimates of the tether properties, 
particularly the damping Cy We have two sources of information as to the 
damping in the tether material, both provided by Martin Marietta Corporation. 
These sources appear to differ substantially. 


Elrst, there is a June 1983 report (Bodley and Park, 1983) which uses a 
value Cy = 9.19 x 10 4 at one point in its analysis (p. 54); the source of this 
number is not clear, and it may be a preliminary value. Second, an internal 
report (Martin Marietta Corporation, 1983) dated September 1983 was kindly 
provided by A. C. Park of MMC; this report gives a "per cent damping" for an 
experimental setup consisting of a mass hanging on the end of a length of 
tether. We have calculated the damping coefficient from the experimental report 
(but see the caveats below) , and converted the earlier number to the same (mks) 


units. In summary: 



Page 33 


- From the experimental data, AE « 10 9 kg-m/s J and C v ~ 200 kg-m/s ("aver- 
aging" the experiments described). These imply a critically damped 
wavelength for tether longitudinal oscillations of about 25 meters (say, 
between 10 and 100) , and a maximum oscillatory frequency of about 30 0 
see -1 . 

- Tre value C v “ 9.19 x 10 4 kg-km/hr given in the "... Orbital Dynamics ..." 
report, when converted to mks system becomes 2.6 x 10 4 kg-m/s. This is 
about 100 times greater than the value calculated from the experimental 
data above. The critically damped wavelength implied is about 3 km and the 
maximum oscillation frequency is about 4 sec* 1 . (The values for AE and p, 
converted to mks, are consistent with the experimental AE and assumed 3 mm 
Kevlar at 1.5 g/cc density.) 

Our interpretation of the experimental results, however, is semawhat uncer- 
tain since (1) the precise testing and computation procedures used were not 
clear, and (2) neither the direct experimental results nor a derived material 
property (C v ) were reported, but a "per cent damping", presumably the damping 
ratio for the experimental setup consisting of a 4.5 kg mass on the end of an 11 
or 21 m "spring." 

In summary, we can say that viscous damping (and possibly other forms of 
energy dissipation such as friction among the strands of the ^arn-like tether) 
can alter the frequencies of longitudinal oscillation modes. The available 
values for the damping coefficient span a wide range, and the maximum oscilla- 
tion frequency corresponding to the larger of the two values is close enough to 
frequencies of interest (i second period or frequency l/2ir) to be of have sub- 
stantial impact. A clearer determination of tether damping properties, par- 
ticularly in a space environment, would be of value. 
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5.0 MISCELLANEOUS EFFORTS 

5.1 Random Vibration Analysis: 

The methods we have been using, and which have been discussed so far in 
this report, are all deterministic in the sense that we consider a given input 
(the atmospheric density as a function of pos^.xon) , make computations with a 
computer simulation program or an analytical model, and observe the result (ac- 
celerations or rotations). The density may be a simple localized enhancement, 
or it may be a complicated randomized function, but during the analysis it is 
fixed; statistical information would need to be garnered from repeated simula- 
tions, or (equivalently under stationarity assumptions) from observing isolated 
sections of a single long simulation. 

Random vibration analysis, which takes into account explicitly the random 
nature of the input, is also likely to be useful. Both the atmospheric density 
variations and the resulting subsatellite acceleration and rotation are random 
processes. The system connecting input and output is still deterministic, but 
rather than simulating the response to an individual realization of the input 
process to achieve an individual realization of the output, we provide a 
statistical description of the input process and gain a statistical description 
of the output process. 

This resulting statistical information complements the specificity of indi- 
vidual case studies. With it we can derive information of use in designing 
instruments, such as the expected time to first excursion to a given response 
level (which might, for instance damage or saturate an instrument) or failure 
due to accumulated stress. Also, the perturbations to the instrumental output 
can then be analyzed to determine their effect on the gravity field after final 
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data analysis. Finally, such considerations might elucidate the knowledge of 
the atmospheric fluctuations required, either from further experiment (or re- 
analysis of old data) or from theoretical arguments (such as gravity wave 
models) . 

Performing such analyses of the dynamically complicated TSS is likely to 
require substantial effort, but the results promise to be of substantial value. 
An early monograph on the subject, lucid but limited, is Crandall and Mark 
(1963) . More complete and recent works are Lin (1967) and Nigam (1983) . 

6.0 DIRECTIONS FOR FUTURE RESEARCH 


As with any major project, the possible tasks seem to increase as work is 
done. We make here some suggestions for future workers in the field. SAO hopes 
to continue dynamic noise studies in connection with the TSS Demonstration Mis- 
sions either under direct contract or in cooperation with other Investigators 
such as Prof. Bergamaschi . 

Improvements in simulation capabilities could proceed along two lines: 

First, a simplified, streamlined and enhanced program similar to SKYHOOK is 
in order. Much ol SKYHOOK'S broac* capabilities are of little importance to 
dynamic noise analysis and only inhibit Investigation of the features of inter- 
est as well as making program maintenance a burden. For situations with rapidly 
varying forcing functions the Gear method may not be the best (other methods 
such as Runge-Kutta-Fehl^erg exist for stiff systems); even if it is, the 
specific subroutine should be tuned to avoid such frequent step size and order 
changes. The system of equations in SKYHOOK is in actuality not merely stiff 
but stiffly oscillatory. There appear to be no general purpose routines or 
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methods for such systems, but there is some promise of developing a special 
purpose method for tethered systems. The data structure should be changed to 
allow use of off-the-shelf band matrix techniques (or special routines written) . 
Subsatellite rotational and aerodynamics should be added, as well as the facil- 
ity for multiple mass subsatellite systems. Rotational effects, particularly 
rotational acceleration, are of key importance in gravity gradiometry. Generat- 
ing equilibrium initial conditions would avoid the current differencing with a 
reference simulation. 

Second, although not a classical form of simulation, the random vibration 
methods should be attempted. The statistical results obtained would provide a 
valuable alternate view to the case studies possible with SKYHOOK like programs. 

For case study simulations, realistic sample atmospheres should be gener- 
ated. Superposition of moderate numbers of local enhancements seems to be of 
limited worth. Methods of generating and using a large grid of (smooth) random 
densities should be explored, including the ad hoc method suggested in Section 
3.3 and the Rotating Lines Method. Also, apart from the question of generating 
random atmospheres, it would be worthwhile to add a damped wave (or set of such) 
to the perturbation options. This feature is suggested by observations of Gross 
and Huang (1964) . 

The tools for analyzing the results of simulation runs can be simplified 
and streamlined. On a more fundamental level, use of the Maximum Entropy Method 
for spectral analysis seems very promising. Both theory and the sxperience 
shown in our figures indicate that we should expect a number of sharp peaks on a 
smoother background. This sort of spectrum is ideally suited to MEM techniques, 
which effectively use rational function approximations to the spectrum rather 
than the simple polynomial approximation of Eourier analysis, and are thus bet- 
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ter able to represent sharp peaks or lines. (Press, et el., 1985.) MEM is 
likely to be useful also for analyzing actual data from a physical gradiometer. 

Design of stabilizing ballast configurations will involve complicated ana- 
lytical studies, probably eventually requiring numerical work for realistic 
designs. Detailed study of mixed partial and ordinary differential systems may 
derive some rules which may be applied to estimate the effects of the continuum 
tether without requiring the full formality. 
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